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Abstract. On a fixed smooth compact Riemann surface with boundary (Mo,g), we show 
that for the Schrodinger operator A + V with potential V 6 C 1,a (Mo) for some a > 0, 
the Dirichlet-to-Neumann map N|r measured on an open set F C dMo determines uniquely 
the potential V. We also discuss briefly the corresponding consequences for potential scat- 
tering at frequency on Riemann surfaces with asymptotically Euclidean or asymptotically 
hyperbolic ends. 



1. Introduction 

The problem of determining the potential in the Schrodinger operator by boundary mea- 
surement goes back to Calderon [8j . Mathematically, it amounts to ask if one can detect some 
data from boundary measurement in a domain (or manifold) with boundary. The typical 
model to have in mind is the Schrodinger operator P := A g + V where g is a metric and V a 
potential, then we define the Cauchy data space by 

C := {(u\ gu ,d u u\gn) ; u G H x (£l), u G kerP} 

where d v is the interior pointing normal vector field to dO,. 

The first natural question is the following full data inverse problem: does the Cauchy data 
space determine uniquely the metric g and/or the potential V? In a sense, the most satisfying 
known results are when the domain Q C M n is already known and g is the Euclidean metric, 
then the recovery of V has been proved in dimension n > 2 by Sylvester-Uhlmann |28| . and 
very recently in dimension 2 by Bukgheim [B] when the domain is simply connected. A re- 
lated question is the conductivity problem which consists in taking V = and replacing A g 
by — divcV where a is a definite positive symmetric tensor. An elementary observation shows 
that the problem of recovering an sufficiently smooth isotropic conductivity (i.e. a = a^Id for 
a function gq) is contained in the problem above of recovering a potential V. For domain of 
M 2 , Nachman [26] used the d techniques to show that the Cauchy data space determines the 
conductivity. Recently a new approach developed by Astala and Paivarinta in [2] improved 
this result to assuming that the conductivity is only a L°° scalar function. This was later 
generalized to L°° anisotropic conductivities by Astala-Lassas-Paivarinta in [3j. We notice 
that there still are rather few results in the direction of recovering the Riemannian manifold 
(£l,g) when V = 0, for instance the surface case by Lassas-Uhlmann [21] (see also [U I15j). 
the real-analytic manifold case by Lassas-Taylor-Uhlmann |20] (see also |14| for the Einstein 
case), the case of manifolds admitting limiting Carleman weights and in a same conformal 
class by Dos Santos Ferreira-Kenig-Salo-Uhlmann [9]. 



i 



2 



COLIN GUILLARMOU AND LEO TZOU 



The second natural, but harder, problem is the partial data inverse problem: if Ti and T% 
are open subsets of <90, does the partial Cauchy data space for P 

e rii r 2 : = {(u\ Vl ,d v u\r 2 ) ; u £ H l (dM Q ), Pu = 0, u = in dtl \ Ti} 

determine the domain f2, the metric, the potential? For a fixed domain of M n , the recovery of 
the potential if n > 2 with partial data measurements was initiated by Bukhgeim-Uhlmann 
[7J and later improved by Kenig-Sjostrand-Uhlmann [18J to the case where T\ and T2 are 
respectively open subsets of the "front" and "back" ends of the domain. We refer the reader 
to the references for a more precise formulation of the problem. In dimension 2, the recent 
works of Imanuvilov-Uhlmann-Yamamoto [16] solves the problem for fixed domains O of M 2 
in the case when Ti = T2 and when the potential are in C 2 > a (Sl) for some a > 0. 

In this work, we address the same question when the background domain is a fixed Riemann 
surface with boundary. We prove the following recovery result: 

Theorem 1.1. Let (Mo,g) be a smooth compact Riemann surface with boundary and let A g 
be its positive Laplacian. For a £ (0, 1), let V±, V2 £ C 1,c *(Mq) be two real potentials and, for 
i = 1,2, let 

(1) S[ =: {(u\ v , cU|r) ; u e H X (M Q ), (A g + Vi)u = 0, u = on dM Q \ T} 

be the respective Cauchy partial data spaces. If Q\ = CJ> then V\ = V2- 

Here the space C^ a {Mo) is the usual ttelder space for a G (0, 1). Notice that when A 9 + Vi 
do not have L 2 eigenvalues for the Dirichlet condition, the statement above can be given in 
terms of Dirichlet-to-Neumann operators. Since A g = e~ 2v A g when g = e 2if> g for some func- 
tion (p, it is clear that in the statement in Theorem ll.lt we only need to fix the conformal class 
of g instead of the metric g (or equivalently to fix the complex structure on Mo). In particu- 
lar, the smoothness assumption of the Riemann surface with boundary is not really essential 
since we can change it conformally to make it smooth and for the Cauchy data space, this just 
has the effect of changing the potential conformally (we only need this new potential to be 
C > a ). Observe also that Theorem 11.11 implies that, for a fixed Riemann surface with boundary 
(Mo,g), the Dirichlet-to-Neumann map on T for the operator u — > —div g ( , yV g u) determines 
the isotropic conductivity 7 if 7 G C 3 ' Q (Mo) in the sense that two conductivities giving rise 
to the same Dirichlet-to-Neumann are equal. This is a standard observation by transforming 
the conductivity problem to a potential problem with potential V := (A g j2)/j2 . So our 
result also extends that of Henkin-Michel [15] in the case of isotropic conductivities. 

The method to reconstruct the potential follows O [TB] and is based on the construction 
of a large set of special complex geometric optic solutions of (A g + V)u = 0, more precisely 
if Tq = 8Mq \ r is the set where we do not know the Dirichlet-to-Neumann operator, we 
construct solutions of the form u = Re(e* ///l (a + r(h)) + e^ e ^ v ^ h s{h) with u|r = where 
h > is a small parameter, $ and a are holomorphic functions on (Mo, g), independent of h, 
H r (MllL 2 = 0(h) while ||s(MIIl 2 = 0(h 3 / 2 \ log h\) as h — > 0. The idea of [6] to reconstruct 
V(p) for p £ Mq is to take $ with a non-degenerate critical point at p and then use stationary 
phase as h — ► 0. In our setting, the function $ needs to be purely real on To an Morse with 
a prescribed critical point at p. One of our main contribution is a geometric construction of 
the holomorphic Carleman weights $ satisfying such conditions. We should point out that 
we use a quite different method than in [16] to construct this weight, and we believe that our 
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method simplifies their construction even in their case. A Carleman estimate on the surface 
for this degenerate weight needs to be proved, and we follows ideas of [16] . We manage to 
improve the regularity of the potential to C 1,a instead of C 2 ' a in [16] . We finally notice 
that we provide a proof, in an Appendix, of the fact that the partial Cauchy data space S r 
determine a potential V E C 0,a (Mo) on T (for a > 0). 

In Section [6l we obtain two inverse scattering results as corollary of Theorem ll.lt first 
for partial data scattering at frequency for A + V on asymptotically hyperbolic surfaces 
with potential decaying at the conformal infinity, and secondly for full data scattering at 
frequenncy for A + V with V compactly supported on an asymptotically Euclidean surface. 

Another straightforward corollary in the asymptotically Euclidean case full data setting is 
the recovery of a compactly supported potential from the scattering operator at a positive 
frequency. The proof is essentially the same as for the operator Aru + V once we know 
Theorem (jl.ip . so we omit it. 

2. Harmonic and Holomorphic Morse Functions on a Riemann Surface 

2.1. Riemann surfaces. We start by recalling few elementary definitions and results about 
Riemann surfaces, see for instance [10J for more details. Let (Mo, go) be a compact connected 
smooth Riemannian surface with boundary dM$. The surface Mq can be considered as a 
subset of a compact Riemannian surface (M,g), for instance by taking the double of Mq 
and extending smoothly the metric go to M. The conformal class of g on the closed surface 
M induces a structure of closed Riemann surface, i.e. a closed surface equipped with a 
complex structure via holomorphic charts z a : U a — > C. The Hodge star operator * acts on 
the cotangent bundle T*M, its eigenvalues are ±i and the respective eigenspace T^qM := 
ker(* + ild) and Tq 1 M := ker(* — ild) are sub-bundle of the complexified cotangent bundle 
CT*M an( j splitting CT*M = T^ M © Tq l M holds as complex vector spaces. Since * is 
conformally invariant on 1-forms on M, the complex structure depends only on the conformal 
class of g. In holomorphic coordinates z = x + iy in a chart U a , one has -k{udx + vdy) = 
—vdx + udy and 



where dz = dx + idy and dz = dx — idy. We define the natural projections induced by the 
splitting of CT*M 



The exterior derivative d defines the De Rham complex — » A — > A 1 — > A 2 — > where A k := 
A k T*M denotes the real bundle of k- forms on M. Let us denote CA k the complexification of 
A fc , then the d and 8 operators can be defined as differential operators d : CA° — ► T^ M and 



T* lfi M\ Ua ~ Cdz, T% tl M\ Ua a Cdz 



7Ti, : CT*M 



Tl M, 7t ,i : CT*M -> T^M. 



d : CA -> T^M by 



(2) 



df := TTifidf, d := n 0t idf, 



they satisfy d = d + d and are expressed in holomorphic coordinates by 



df = d z f dz, df = dzf dz. 
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with d z := \(d x — id y ) and <9 ? := \{d x + id y ). Similarly, one can define the d and d operators 
from CA 1 to CA 2 by setting 

d(u 1}0 + o; ,i) := dwo,i, d(ui >0 + u ,i) := du lfi 
if cjo,i G Tq l M and u)\ o G T{* M. In coordinates this is simply 

d{udz + udz) = dv A dz, d(udz + wdz) = du A dz. 
There is a natural operator, the Laplacian acting on functions and defined by 

A/ := -2i*8df = d*d 

where d* is the adjoint of d through the metric g and * is the Hodge star operator mapping 
A 2 to A and induced by g as well. 

2.2. Maslov Index and Boundary value problem for the d Operator. In this subsec- 
tion we consider the setting where M is an oriented Riemann surface with boundary OMq and 
r C 8Mq is an open subset and we let Tq = OMq \ V be its complement in 8Mq. Following 
[22], we adopt the following notations: let E — » Mq be a complex line bundle with complex 
structure J : E -» E and let D : C°°(M ,£) -> C°°(Mo, Tq X ® £) be a Cauchy-Riemann 
operator with smooth coefficients on Mo, acting on sections of the bundle E. Observe that 
in the case when E = Mq x C is the trivial line bundle with the natural complex structure 
on Mo, then D can be taken as the operator d introduced in ([2]). For q > 1, we define 

D F : W e /(M ,E) - W e ~ 1 ' q (M ,TQ 1 M (8) E) 

where F C E \qm is a totally real subbundle (i.e. a subbundle such that JF n F is the zero 
section) and Dp is the restriction of D to the L 9 -based Sobolev space will £ derivatives and 
boundary condition F 

W e /(M ,E) := U G VF^(M ,£) | e(9M ) C F}. 
In this setting, we have the following boundary value Riemann- Roch theorem stated in |22| : 

Theorem 2.1. Let E — > Mq be a complex line bundle over an oriented compact Riemann 
surface with boundary and F C E \dM be a totally real subbundle. Let D be a smooth Cauchy- 
Riemann operator on E acting on W e,q (Mo, E) for some q > 1 and t G N. Then 

1 ) The following operators are Fredholm 

Dp : W £ /(M , E) - W'-^Mq, T^M ® F) 

F> : W% q (M ,lZ tl M E) W'-^iM^E). 

2) The real Fredholm index of Dp is given by 

lnd(D F ) = x(M ) + ^E,F) 

where x(Mo) is the Euler characteristic of Mq and /u(F, F) is the boundary Maslov index of 
the subbundle F . 

3) If n(E,F) < 0, then Dp is injective, while if fi(E,F) + 2%(Mo) > the operator Dp is 
surjective. 
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In the case of a trivial bundle E = Mq x C and 8Mq = UJL^Mq is a disjoint union of 
m circles, the Maslov index can be defined (see [22j p. 554-555]) to be the degree of the map 
p o A : 8M -» BM where 

AkM : S 1 ~ 8 l M -> GL(1,C)/GL(1,R) 

is the natural map given by the totally real subbundle (the space GL(1,C)/GL(1,R) being 
the space of totally real subbundles of C over a circle) and 

p : GL(1,C)/GL(1,M) -► S 1 , p(AGL(l,R)) := det(A 2 )/det(AM). 

As an application, we obtain the following (here and in what follows, H m (Mo) := W m,2 {MQ)): 

Corollary 2.2. (i) For q > 1 and k G N , let uj G W k ' q (M , Tq^Mq), then there exists 
u G iy fc+1,,3, (Mo) holomorphic on Mq, real valued on Tq, such that Bu = u. 
(ii) Form > 1/2, let f G H m (dMo) be a real valued function, then there exists a holomorphic 
function v G i/ m+ 2(M ) such that Re(w)|r = /■ 

(Hi) For k G N and g > 1, i/te space of W k ' q (Mo) holomorphic functions on Mq which are 
real valued on Tq is infinite dimensional. 

Proof, (i) Let L G N be arbitrary large and let us identify the boundary as a disjoint union 
of circles BMq = WJL i diM$ where each 9jMo ~ S 1 . Since F will be the piece of the bound- 
ary where we know the Cauchy data space, it is sufficient to assume that T is a connected 
non-empty open segment of 8\Mq = S , and which can thus be defined in a coordinate 6 
(respecting the orientation of the boundary) by T = {6 G S* 1 | < 8 < 2ir/k} for some integer 
k. Define the totally real subbundle of F C E\qm = ]J l j=i(9jMo x C) by the following: on 

d 1 M ~ S 1 parametrized by 9 G [0, 2vr], define F e = e ia(e) R. Where a : [0, 2vr] -> R is a smooth 
nondecreasing function such that a(9) = in a neighbourhood [0,e] of 0, a(2ir/k) = 2Lir for 
some L G N, and a{9) = 2Lir for all 9 > 2tt jk. For the rest of diM§, ..,d m Mo, we just let 
F\diMa = 5" 1 x R. The Maslov index p(E,F) is then given by 2L and so, by theorem 12. 11 Dp 
is surjective if 2x(Mq) + 2L > 0. Since L can be taken as large as we want this achieves the 
proof of (i). 

(ii) Let w G H m+ 2 (Mq) be a real function with boundary value / on 8Mq, then by (i) there 
exists R G -fT m+1 / 2 (Mo) such that idR = —Bw and R purely real on To, thus v := iR + w is 
holomorphic such that Re(u) = / on Tq. 

(iii) Taking the subbundle F as in the proof of (i), we have that dimker Dp = x(Mq) + 2L 
if L satisfies 2x(Mq) + 2L > 0, and since L can be taken as large as we like-, this concludes 
the proof. □ 

Lemma 2.3. Let {po,p±, .-,Pn} G Mq be a set of n + 1 disjoint points. Let c\, . . . ,ck G C, 
N G N, and let z be a complex coordinate near pq such that pq = {z = 0}. Then if po G 
int(Mo), there exists a holomorphic function f on Mq with zeros of order at least N at each 
Pj, such that f is real on Tq and f(z) = cq + c\Z + ... + ckz k + 0(|z| ft:+1 ) in the coordinate 
z. If pq G BMq, the same is true except that f is not necessarily real on Tq. 

Proof. First, using linear combinations and induction on K, it suffices to prove the Lemma 
for any K and cq = ■ ■ ■ = ck-\ = 0, which we now show. Consider the subbundle F as in the 
proof of (i) in Corollary 12.21 The Maslov index p(E, F) is given by 2L and so for each G N, 
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one can take L large enough to have fi(F,E) + 2x(Mq) > 2N(1 + n). Therefore by Theorem 
12,11 the dimension of the kernel of Of will be greater than 2(n + l)iV. Now, since for each 
Pj and complex coordinate Zj near pj, the map u — > (u(pj),d Zj u(pj), . . . , d^~ l u{pj)) G is 
linear, this implies that there exists a non-zero element u £ keiDp which has zeros of order 
at least iV at all pj . 

First, assume that pq £ int(Mo) and we want the desired Taylor expansion at pq in the 
coordinate z. In the coordinate z, one has u(z) = az M + 0(\z\ M+1 ) for some a / and 
M > N. Define the function tk{z) = xi z )^~ z ~ M+K where x{ z ) is a smooth cut-off function 
supported near p$ and which is 1 near pq = {z = 0}. Since M > N > 1, this function 
has a pole at po and trivially extends smoothly to Mq\{pq}, which we still call tk- Observe 
that the function is holomorphic in a neighbourhood of pq but not at pq where it is only 
meromorphic, so that in Mq \ {po}, drx is a smooth and compactly supported section of 
Tq 1 Mq and therefore trivially extends smoothly to Mq (by setting its value to be at po) to 
a one form denoted u>k- By the surjectivity assertion in Corollary 12,21 there exists a smooth 
function Rk satisfying dRx = — and that Rk\f 6 ^- We now have that Rk + rx is 
a holomorphic function on M\{po} meromorphic with a pole of order M — K at po, and in 
coordinate z one has z M ~ K (Rk(z) + rxiz)) = ck + 0(|z|). Setting fx = u{Rk + rj<), we 
have the desired holomorphic function. Note that / also vanish to order TV at all p%, . . . ,p n 
since u does. This achieves the proof. 

Now, if pq £ OMq, we can consider a slightly larger smooth domain of M containing Mq 
and we apply the the result above. □ 



2.3. Morse holomorphic functions with prescribed critical points. The main result 
of this section is the following 

Proposition 2.1. Let p be an interior point of Mq and e > small. Then there exists a 
holomorphic function $ on Mq which is Morse on Mq (up to the boundary) and real valued on 
Tq, which has a critical point p' at distance less than e from p and such that Im(<&(p')) ^ 0. 

Let be a connected open set of M such that is a smooth surface with boundary, 
with Mq C and Tq C 80. Fix k > 2 a large integer, we denote by C fc (0) the Banach 
space of C k real valued functions on 0. Then the set of harmonic functions on which 
are in the Banach space C k (0) (and smooth in by elliptic regularity) is the kernel of the 
continuous map A : C k (0) — » C fc ~ 2 (0), and so it is a Banach subspace of C k (0). The set 
IK C C k (0) of harmonic functions u in C k (Q) such there exists v S C k (0) harmonic with 
u + iv holomorphic on is a Banach subspace of C k {0) of finite codimension. Indeed, let 
{71, ..,7tv} be a homology basis for 0, then 

"K = kerL, with L : ker A n C k (0) -> defined by L(u) := (J^ J duj ^ ^ 

For all r' C 8Mq such that the complement of T' contains an open subset, we define 

K n := {uSM;u| n =0}. 

We now show 

Lemma 2.4. The set of functions u G ^r' which are Morse in is residual (i.e. a countable 
intersection of open dense sets) in "K T t with respect to the C k (0) topology. 
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Proof. We use an argument very similar to those used by Uhlenbeck [30]. We start by 
denning m : x iHp^ — > T*0 by (p,u) i— ► (p,du(p)) G TpO. This is clearly a smooth map, 
linear in the second variable, moreover m u := m(., u) = (•, du{-)) is Fredholm since is finite 
dimensional. The map u is a Morse function if and only if m u is transverse to the zero section, 
denoted T *O, of T*0, ie. if 

(T *O) = T mti(p) (T*0), Vp G such that m«(p) = (p, 0), 

which is equivalent to the fact that the Hessian of it at critical points is non-degenerate (see 
for instance Lemma 2.8 of [3D]). We recall the following transversality theorem ( \S0\ Th.2]): 

Theorem 2.5. Let m : X x 3f r / —> W be a C k map, where X, "K^ , and W are separable 
Banach manifolds with W and X of finite dimension. Let W' C W be a submanifold such 
that k > max(l,dimX — dimPF + dim W'). If m is transverse to W then the set {u G 
JCp^ ; m u is transverse to W'} is dense in JC^ , more precisely it is a residual set. 

We want to apply it with X := 0, W := T*0 and W 1 := Tq 0, and the map m is defined 
above. We have thus proved Lemma 12.41 if one can show that m is transverse to W . Let 
0,it) such that m(p,u) = (p,0) G W. Then identifying T (p0) (T*0) with T p 0®T*0, one has 

D(p U \m(z,v) = (z,dv(p) +Hess p (n)z) 

where Hess p u is the Hessian of u at the point p, viewed as a linear map from T p to T*0. 
To prove that m is transverse to W' we need to show that (z,v) — > (z,dv(p) + H.ess p (u)z) is 
onto from T p © ^r' to T p (BT*0, which is realized for instance if the map v — > dv{p) from 
J{p' to T*0 is onto. But from Lemma 12.31 we know that there exist holomorphic functions 
v and v on such that v and v are purely real on T' . Clearly the imaginary parts of v and 
v belong to Hr' - Furthermore, for a given complex coordinate z near p = {z = 0}, we can 
arrange them to have series expansion v(z) = z + 0(|z| 2 ) and v(z) = iz + 0(|z| 2 ) around the 
point p. We see, by coordinate computation of the exterior derivative of Im(u) and Im(-u), 
that dlm(v)(p) and dlm(v)(p) are linearly independent at the point p. This shows our claim 
and ends the proof of Lemma 12.41 by using Theorem 12.51 □ 

We now proceed to show that the set of all functions u G Jtjv such that it has no degenerate 
critical points on T' is also residual. 

Lemma 2.6. For all p G T' Q and k G N, there exists a holomorphic function u G C k (0), such 
that Im(it)| r ^ = and du(p) / 0. 

Proof. The proof is quite similar to that of Lemma 12.31 By Lemma 12.31 we can choose a 
holomorphic function v G C fc (0) such that v(p) = and Im(i;)|p/ = 0, then either dv(p) ^ 
and we are done, or dv(p) = 0. Assume now the second case and let M G N be the order 
of p as a zero of v. By Riemann mapping theorem, there is a conformal mapping from 
a neighbourhood U p of p in to a neighbourhood {\z\ < e,Im(z) > 0} of the real line 
Im(z) = in C, and one can assume that p = {z = 0} in these complex coordinates. 
Take r(z) = x{z)z~ M+1 where \ G Co°(|z| < e) is a real valued function with x( z ) = 1 
in {\z\ < e/2}. Then Br vanishes in the pointed disc < \z\ < e/2 and it is a compactly 
supported smooth section of Tj" O outside, it can thus be extended trivially to a smooth 
section of T^ O denoted by u. We can then use (i) of Corollary 12.21 there is a function R 
such that dR = —to and Im(i?)|jv = 0, and so B{R + r) = in \ {p} and R + r is real valued 
on Tq (remark that r is real valued on T' ) and has a pole at p of order exactly M — 1. We 
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conclude that u := v(R + r) satisfies the desired properties, it vanishes at p but with non zero 
complex derivative at p. □ 

Lemma 2.7. Let T' Q C 80 be an open set of the boundary. Let <fi : — > R be a harmonic 
function with 4>\^i = 0. Let p £ T' Q be a critical point of <j), then it is nondegenerate if and 
only if d T d u u ^ where d T and d u denote respectively the tangential and normal derivatives 
along the boundary. 

Proof. By Riemman mapping theorem, there is a conformal transformation mapping a 
neighbourhood of p in to a half-disc D := {\z\ < e,Im(z) > 0} and <90 = {Im(z) = 0} 
near p. Denoting z = x + iy, one has (c^ + dy)<p = in D and dx4>\ y =o = 0, which implies 
dy(j)(p) = 0. Since d v = d y and d T = d x for some smooth function /, and since d<j){p) = 0, 
the conclusion is then straightforward. □ 

Let N*dO be the conormal-bundle of 80 and A r *TQ be the restriction of this bundle to T' . 
Denote the zero sections of these bundles respectively by N^dO and NqT'q. We now define 
the map 

b : r' X ^ -> A^*^, b{p, u) := (p, d u u). 
For a fixed u G ^r' : we a l so define b u {-) := b(-,u). Simple computations yield the 

Lemma 2.8. Suppose that p G T' Q is such that d u u(p) = 0, then d T d v u{p) ^ if and only if 

Imaged) + ^W^Io) = T {pfi) (N*T' ). 
Proof. This can be seen by the fact that for all p £ T' Q such that b u {p) = (p, 0), 

D p b u : T p T' -» T {Pi0) (N*T' ) ~ T p V N;V 
is given by w i— > (w,d T d u u(p)w). □ 



At a point (p, tt) such that u) = 0, a simple computation yields that the differential 
D ( P ,u)b ■ T p T' x IK r / o -> r (Pi ^ n{p)) (iV*r / ) is given by ^ (w,d T d u u(p)w + d v u'(jp)). 

This observation combined with Lemma 12.61 shows that for all (p, u) £ T'q x IKp' such that 
b(p,u) = (p, 0), 6 is transverse to ATqTq at (p, 0). Now we can apply Theorem 12.51 with 
X = T' Q , W = N*T' Q and W = N$T' we see that the set {u £ "K v ,;b u is transverse to 
is residual in !Kjv • In view of Lemmas 12.7} we deduce the 

Lemma 2.9. The set of functions u £ "K^/ such that u has no degenerate critical point on 
T'q is residual in "K^t . 

Observing the general fact that finite intersection of residual sets remains residual, the 
combination of Lemma 12.91 and Lemma 12.41 yields 

Corollary 2.10. The set of functions u £ ^r' which are Morse in and have no degenerate 
critical points on T' is residual in Jir' with respect to the C k {0) topology. In particular, it 
is dense. 

We are now in a position to give a proof of the main proposition of this section. 

Proof of Proposition 12.11 As explained above, choose in such a way that is a smooth 
surface with boundary, containing Mq, that Tq C 80 and contains dMo\To. Let T' be 
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an open subset of the boundary of such that the closure of Tq is contained in T' Q and 
<90\rQ ^ 0. Let p be an interior point of Mq. By lemma [2T3| there exists a holomorphic 
function / = u + iv on such that / is purely real on r' , v(p) = 1, and df(p) = (thus 

V G Ji T r). 

1 o' 

By Corollary I2.1UI there exist a sequence {vj)j of Morse functions Vj G ^r' such that 
Vj — > f in C fc (Mo) for any fixed A; large. By Cauchy integral formula, there exist harmonic 
conjugates Uj of Vj such that := Uj + ifj — > / in C fc (Mo). Let e > be small and let 
U C be a neighbourhood containing p and no other critical points of f, and with boundary 
a smooth circle of radius e. In complex local coordinates near p, we can identify df and dfj 
to holomorphic functions on an open set of C. Then by Rouche's theorem, it is clear that dfj 
has precisely one zero in U and Vj never vanishes in U if j is large enough. 

Fix $ to be one of the fj for j large enough. By construction, <I> is Morse in and has 
no degenerate critical points on Tq C T' q . We notice that, since the imaginary part of $ 
vanishes on all of T' , it is clear from the reflection principle applied after using the Riemann 
mapping theorem (as in the proof of Lemma 12. 7p that no point on Tq C T' can be an accu- 
mulation point for critical points. Now 8Mq\Tq is contained in the interior of and therefore 
no points on <9Mo\ro can be an accumulation point of critical points. Since $ is Morse in 
the interior of 0, there are no degenerate critical points on 9Mo\rQ. This ends the proof. □ 



3. Carleman Estimate for Harmonic Weights with Critical Points 

In this section, we prove a Carleman estimate using harmonic weight with non-degenerate 
critical points, in way similar to [16] . Let us define Tq := {p G 8Mq; d u ip(p) = 0} and let 
T := 8Mq \ Tq its complement. 

Proposition 3.1. Let (MQ,g) be a smooth Riemann surface with boundary, and let ip : Mq — > 

M. be a C k (Mo) harmonic Morse function for k large. Then for all V G L°°(Mq) there exists 
an }iq > such that for all h G (0, ho) and u G C°°(M) with u\qm = 0, we have 

1 1 

^IMIl2(M ) + -tfW U \ dl P\\\h(Mo) + W du \\h(M ) + W d » u \\h(r ) 

< C(\\e-^(A 9 + V)e^ h u\\l HMo) + ^IIMlS^) 
where d v is the exterior unit normal vector field to 8Mq. 

Proof. We start by modifying the weight as follows: if ipQ := ip : Mq — > M. is a real val- 
ued harmonic Morse function with critical points {p±, . . . ,pn} in the interior of Mq, we let 
ip* : Mq — > R be harmonic functions such that pj is not a critical point of <pj for j = 1, . . . , N , 
their existence is insured by Lemma 12.31 For all e > 0, we define the convexified weight 

^:=^-£(£f=oN 2 )- 

To prove the estimate, we shall localize in charts £lj covering the surfaces. These charts 
will be taken so that if Qj n BMq ^ 0, then Qj n 8Mq ~ S 1 is a connected component of 8Mq. 
Moreover, by Riemann mapping theorem (e.g. Lemma 3.2 of [23]), this chart can be taken to 
be a neighbourhood of \z\ = 1 in {z G C; \z\ < 1} and such that the metric g is conformal to 
the Euclidean metric \dz\ 2 . 
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Lemma 3.1. Let VL be a chart of Mq as above and (p e : £1 ^ M. be as above. Then there are 
constants C,C' > such that for all u G C°°(M) supported in £1 and h > small enough, 
the following estimate holds: 
(4) 

-Mi^o) + C '{ ~ M(drU,u) L 2 idMo) ) + \ [ \u>\ 2 d v <p t dv g ) < \\e-^ h de^ h u\\l, [Mo) 



€ 



where d u and d T denote respectively the exterior pointing normal vector fields and its rotation 
by an angle +n/2. 

Proof. We use complex coordinates z = x + iy in the chart O where u is supported. Observe 
that the Lebesgue measure dxdy is bounded below and above by dv 9 , g is conformal to \dz\ 2 
and the boundary terms in (JH) depend only on the conformal class, it suffices to prove the 
estimates with respect to dxdy and the Euclidean metric. We thus integrate by parts with 
respect to dxdy 



A\\e-^l h de^ h uf =|| + i ^)u + (id v + ^)u\ 

h y h 



2 



Ufa , id V^ \ ||2 . , \ l|2 

=\\{d x + —^—)u\\ + \\{idy + —f^Ml 

+ \ j ^feH 2 -^d x <f e .d x \u\ 2 -^dyif e .d y \u 

(5) 2 f f / 

+ — / c^eM 2 — 2 / ( <9 x Re(u).<9 y Im(u) — 9 a .Im(u).5j / Re(ti) 

= 11(9, + ^H| 2 + \\{id y + ^)uf + IJ^ A^\u\ 2 

+ T [ d v tp e \u\ 2 + 2 / (9 T Re(u).Im(n). 

where A := —{d 2 + d 2 ), d u is the exterior pointing normal vector field to the boundary and 
d T is the tangent vector field to the boundary (i.e. d v rotated with an angle vr/2) for the 
Euclidean metric \dz\ 2 . Then (uAip e ,u) = ^(\d(fo\ 2 + \d(pi\ 2 + .. + |c^jv| 2 )M 2 , since fj are 
harmonic, so the proof follows from the fact that |d</9 | 2 + (d^il 2 + •• + I^aH 2 is uniformly 
bounded away from zero. □ 



The main step to go from @ to © is the following lemma which is a slight modification 
of the proof in [TB] : 

Lemma 3.2. With the same assumptions as in Proposition \3.1[ and ifVt is either an interior 
chart of (M, g) or a chart containing a whole boundary connected component, then there are 
positive constants c and C such that for all e > small, all < h -C e and all smooth function 
u supported in O with u\qm = 0, we have 

c(||e-^/%e^||| 2(A/o) + \\\d„u\\h { r)) > 

c / 1 1 1 \ 

-(l\\ U \\h(M) + tf\\ u \ d( P\\\l*{M) + tfhldfMhiM) + \\ du \\h(M)J + \\ d Ml*(dM ) 

where d v denote the unit normal vector field to 8Mq. 
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Proof. Since the norms induced by the metric g in the chart are conformal to Euclidean 
norms, and there is / smooth such that A g = — e 2 ^ (d 2 +d 2 ) = e 2 $ A in the complex coordinate 
chart, it suffices to get the estimate © for Euclidean norms and Laplacian. Clearly, we can 
assume u G Hq(M) to be real valued without loss of generality. Now let Q(z) be a holomorphic 
function in f2 which has no zeros, then by (0) we have 

\\Q e -^ /h Ae^ /h u\\ 2 = l§\\e-^ lh de^ /h Qe-^ /h de^ /h u\\ 2 > 

— \\du + ^ull 2 - 4Im( / drUJ.uo) + \ f \u\ 2 d v <p e 
eh V J dMo ) h J 9Mo 

with u := Qe~ lfic ^ h de !fit ^ h u and here everything is measured with respect to Euclidean metric 
dx 2 + dy 2 and measure dxdy. Since u\qm = 0, one has uj\qm = {A + iB)d u u where A + iB = 
Q((d u ,d x ) — i(d u ,d y )) and — hn(d T uj.u)) = (d T A.B — Ad T B)\d u u\ 2 we deduce that for some 
c> 

\\Q e -^/h Ae ^/h u \\2 > 



( 7 ) -(J|du|| + ^ ||it|dyj e ||| + ~(d x u,ud x <p e ) + -(dyu,udy(p. 

+4 / (Ad T B-d T AB)\d u u\ 2 + t I \Q\ 2 \d u u\ 2 d»Ve- 

JdM h JdM 

Using the fact that u is real valued, that ip is harmonic and that ^2f = o \ d<Pj\ 2 is uniformly 
bounded away from 0, we see that 

^(d x u,ud x ip e ) + ^(d y u,udyip t ) = ~{u,uA(p e ) > ^||u|| 2 
for some C > and therefore, 

(8) \\Qe-^ /h Ae^ /h u\\ 2 > -(\\du\\ 2 + ^\\u\dp e \\\ 2 + -\\u\\ 2 ) + boundary terms. 

e h z e 

Now if the diameter of the support of u is chosen small (with size depending only on 
|Hessy?o|(p)) with a unique critical point p of ipo inside, one can use integration by parts 
and the fact that the critical point is non-degenerate to obtain 

(9) \\du\\ 2 + ^\\u\d^\\\ 2 >l 



dz(u )d z ip dxdy 



1 

> - 
- h 



u d 2 ifo dxdy 



> T NI 



for some C > 0. Clearly the same estimate holds trivially if f2 does not contain critical point 
of tpo. Using a partition of unity (9j)j in Q and absorbing terms of the form ||u<9#j|| 2 into 
the right hand side, one obtains Q for any function u supported in and vanishing at the 
boundary. Thus, combining with (jSJ), there are positive constants c,c',C" such that for h 
small enough 

^(||dn|| 2 + ^\\ u \dcp e \\\ 2 + j\\u\\ 2 ) > - e (\\du\\ 2 + ^||u|^ ||| 2 - ^\\u\\ 2 ) 

> -(\\duf + ^\\u\d<p \\\ 2 + ~\\u\\ 2 ). 
Combining now with (JS]) and using that \Q\ is bounded below gives 

\\ e -<p a /h Ae <p./h u u2 > _(Mu|| 2 + -L||n|^||| 2 + rlkll 2 ) + boundary terms. 
e h z h 



12 COLIN GUILLARMOU AND LEO TZOU 

Let us now discuss the boundary terms in ([?]). If ipj are taken so that d v (fj = on Tq, then 
d u ip e = on To and d u (p e = d u ip + 0(h/e) on T and thus 

\ [ \Q\ 2 \d v u\ 2 \d v ipe\ < ^ / \d„u\ 2 

n JdMo 11 JV 

for some constant C\. We finally claim that there exist Q with no zeros in $7 such that 
Ad T B — Bd T A is bounded below by a positive constant on 8Mq n fi. Indeed, since the 
chart near a connected component can be taken to be an interior neighbourhood of the circle 
\z\ = 1 in C, one can take A + iB = e J * where i 6 S 1 parametrize the boundary compo- 
nent, so that Ad T B — Bd T A = 1 since d T = dt for the Euclidean metric. Since moreover 
(d v , d x ) — i(d v , d y ) = A — iB = e~ lt , we deduce that on the boundary Q(t) = e 2tt and so it 
suffices to take Q(z) = z 2 . This achieves the proof. □ 

Proof of Proposition [3TT1 Using triangular inequality and absorbing the term ||yu|| 2 into 
the left hand side of ((3j) , it suffices to prove ([3]) with A g instead of A g + V. Let v G Cg°(M), 
we have by Lemma 13.21 that there exist constants c, d ,C,C > such that 

K^" e ~ 1P ' /ht '" 2 + ^ e ~^ /kv \ d v\W 2 + ■jp\\ e ~' P ' ,hv \<hPe\f + He-^dwll 2 ) + \\<r* e/h d u v\\l 



< 

3 



E^^II^Vll 2 + ^\\e-^ lh XjvW\\\ 2 + ±\\ e -^ h Xl v\dip t 
+ l|e-^ / Mx^)H 2 )+l|e-^ / ^^||r 



< C( £ \\e~^/ h A g ( Xj v)\\ 2 + \\e-^ h d u v\\ 2 r 



3 



< C'(\\e-^ h A g v\\ 2 + \\e-^ /h v\\ 2 + \\ e -^ h dvf + \\e-^ /h d u v\\l) 

where (xj)j is a partition of unity associated to the complex charts 0,j on M. Since constants 
on both sides are independent of e and h, we can take e small enough so that C'\\e-^/ h v\\ 2 + 
C'\\e~' P€ / h dv\\ 2 can be absorbed to the left side. Now set v = e Ve f h w with w\qm = 0) then 
we have 

-||U>|| 2 + TjjHH^I IP + ■pll u 'l^e||| 2 + \\dw\\ 2 + ll^ w llr 

< c{\\e^l h A g e^l h w\\ 2 + \\d v u\\i) 

Finally, fix e > and set u := e e ^j^ w and use the fact that e e is independent 

of h and bounded uniformly away from zero and above, we then obtain the desired estimate 
for < h <C e. □ 



4. Complex Geometric Optics on a Riemann Surface 

As explained in the Introduction, the method for identifying the potential at a point p is 
to construct complex geometric optic solutions depending on a small parameter h > 0, with 
phase a Carleman weight (here a Morse holomorphic function), and such that the phase has 
a non-degenerate critical point at p, in order to apply the stationary phase method. In this 
section, the potential V has the regularity V E C ,a (Mo) for some a > 0. 
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Choose p G int(Mo) such that there exists a holomorphic function $ = cp + itp which is 
Morse on Mq, C k in Mq for large k G N and such that d&(p) = and <3? has only finitely 
many critical points in Mq. Furthermore we ask that $ is purely real on Tq. By Proposition 
12, II such points p form a dense subset of Mq. Given such a holomorphic function, the purpose 
of this section is to construct solutions u on Mq of (A + V)u = of the form 

(10) u = e q>/h (a + ha + r 1 ) + e*/ h (o + ha + n) + e f/h r 2 with ii|r = 

for h > small, where a is holomorphic and u G C fc (Mo) for large G N, ao G H 2 (Mq) 
is holomorphic, moreover a(p) 7^ and a vanishes to high order at all other critical points 
p' G Mq of Furthermore, we ask that the holomorphic function a is purely imaginary on 
Tq. The existence of such a holomorphic function is a consequence of Lemma 12.31 Given 
such a holomorphic function on Mq we consider a compactly supported extension to M, still 
denoted a. 

The remainder terms n, T2 will be controlled as /i — > and have particular properties near 
the critical points of $. More precisely, r2 will be a 0^2 (/i 3 / 2 | log h\) and n will be of the 
form hr\2 + °l 2 (M where fxi is independent of h, which can be used to obtain sufficient 
informations from the stationary phase method in the identification process. 

4.1. Construction of r\. We shall construct r\ to satisfy 

e~^ h (A g + V)e^ h (a + n) = L z(h\ ]ogh\) 

and ri = m + /iri2- We let G be the Green operator of the Laplacian on the smooth surface 
with boundary Mq with Dirichlet condition, so that A 9 G = Id on L 2 (Mq). In particular this 
implies that BdG = where * -1 is the inverse of * mapping functions to 2-forms. We 

extend a to be a compactly supported C k function on Mq and we will search for r\ G H 2 (Mq) 
satisfying ||ri|| L 2 = 0(h) and 

(11) e -^/h de 2^/h n = -dG(aV) + uj + H i(h\\ogh\) 

where a; is a smooth holomorphic 1-form on Mq. Indeed, using the fact that $ is holomorphic 
we have 

e-*/ h A g e*/ h = -U*de-*l h de*l h = -2i * Be'^-^de^-^ = -2i * Be'^de 2 ^ 

and applying —li-kd to (fTTI) . we obtain (note that dG(aV) G C 2,Q (Mo) by elliptic regularity) 

e"*/ h (A fl + V>*/ fc n = -aV + L2 (h\ log h\). 

We will choose w to be a smooth holomorphic 1-form on Mq such that at all critical point p 1 
of <!> in Mq, the form b := <9G(aF) — u with value in T^qMq vanish to the highest possible 
order. Writing b = b(z)dz in local complex coordinates, b(z) is C 2+a by elliptic regularity 
and we have —2idzb(z) = aV, therefore d z dgb(p') = d 2 b(p') = at each critical point p' 7^ p 
by construction of the function a. Therefore, we deduce that at each critical point p' ^ p, 
8G(aV) has Taylor series expansion Ylj=o c j z ^ + 0(|2| 2+a ) for some N large. That is, all the 
lower order terms of the Taylor expansion of dG(aV) around p' are polynomials of z only. 

Lemma 4.1. Let {po, ~.,pn} be finitely many points on Mq and let 6 be a C 2,a section of 
T^qMq. Then there exists a C k holomorphic function f on M with k G N large, such that 
u = df satisfies the following: in complex local coordinates z near pj , one has d l z 6(pj) = 
d^uj(pj) for £ = 0, 1, 2, where 9 = 8(z)dz and uj = oj(z)dz. 
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Proof. This is a direct consequence of Lemma 12.31 □ 

Applying this to the form dG{aV) and using the observation we made above, we can construct 
a C k holomorphic form oj such that in local coordinates z centered at a critical point p' of 
(i.e p' = {z = 0} in this coordinate), we have for b = dG(aV) — to = b(z)dz 

\dfdib{z)\=0{\z\ 2+a - t - m ), for^ + m<2, ifp'/p 
\b(z)\=0(\z\), iip>=p. 

Now, we let xi S C^°(Mq) be a cutoff function supported in a small neighbourhood U p of 
the critical point p and identically 1 near p, and \ G C^°(Mq) is defined similarly with \ = 1 
on the support of \i- We will construct r\ = rn + hr\ 2 in two steps : first, we will construct 
Tix t° solve equation (JTTj) locally near the critical point p of and then we will construct the 
global correction term r\ 2 away from p by using the extra vanishing of b in (|12|) at the other 
critical points. 

We define locally in complex coordinates centered at p and containing the support of % 

r n ■= xe- 2 ^ lh R{e 2l1p/h xib) 

where Rf{z) := — (27ri) _1 J R2 jz^fd£ A d£ for / G L°° compactly supported is the classical 
Cauchy-Riemann operator inverting locally d z (rn is extended by outside the neighbourhood 
of p). The function rn is in C 3+a (Mo) and we have 



-2* 



iV/^ (e 2#/V u ) = Xl (-dG(aV) + u) + r, 



(13) 

with v := e- 2i ^ lh R{e 2 ^ /h xib)d X . 

We then construct r\i by observing that b vanishes to order 2 + a at critical points of $ other 
than p (from (I12p ). and d\ = in a neighbourhood of any critical point of ip, so we can find 
r\2 satisfying 

2ir 12 dil) = (1 - Xi)b- 

This is possible since both dip and the right hand side are valued in Tj* Mo, dip has finitely 
many isolated on Mq: t\i is then a function which is in C 2,01 (Mq\P) where P := {p±, . . . ,pn} 
is the set of critical points other than p, it extends to a C 1,a (Mo) and it satisfies in local 
complex coordinates z near each pj 

\dldjr 12 {z)\ <C\z-p j \ 1+a - /3 -~<, /3 + 7 <2. 

by using also the fact that dip can be locally be considered as holomorphic function with a 
zero of order 1 at each pj. This implies that r\ G H 2 (Mq) and we have 

e -m/h Q ^/h r ^ = b + hdri2 = _ 5G ( aF ) _ u + hdri2 + r? . 

Now the first error term | \dri2\ \h 1 (m ) i s bounded by 

:i-xiw*)i 



\\dri 2 \\m(M ) < C 



h. To deal with the rj term, we need the following 



d z ip(z 

i_ 



< c 



H 2 (U p ) y 



for some constant C, where we used the fact that ^g^r^ is in H 2 (U P ) and independent of 
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Lemma 4.2. The following estimates hold true 
|Mlfl* = 0(|log/i|), \\ V \\ H i<0(h\logh\), 
where ryi solves Ivryidty = b. 

Proof. We start by observing that 



|n|| L 2 = 0(h), \\r 1 - hr 



12\\L 2 



o(h) 



(14) 



-2iip/h jjf 2iip/h 



R{e zt ^ n Xl b) ~ h 



L\U P ) 



+ h\\r 12 \ 



\n - hr 12 \\ L 2 



< 



Xe 



-2H,/h R ( e w/h Xl Q _ h 



Xib 



28^ 



L 2 (M ) + h 

T](Z 



rj(z) 



2d z ^ 



L 2 (U P ) 



+ h 



L 2 (U P ) 



2d z ^ 



L 2 (U P ) 



and we will show that the first and last terms in the right hand sides are o(h) while h\ |n.2| |l 2 (m ) 
is 0(h). The first term is estimated in Proposition 2.7 of [16], it is a o(h), while the last term 
is clearly bounded by C/i||?7||^2 and the middle one by Ch for some constant C by using that 
dijj does not vanish on the support of r] and the fact that b vanishes at critical points of ip. 
Now are going to estimate the H 2 norms of r\. Locally in complex coordinates z centered at 
p (ie. p = {z = 0}), we have 



(15) 



rj(z) 



-d z x(z)e h 



e h ■ 



7T 



£ = £i + 



Since b is C 2,a in U, we decompose 6(£) = (V6(0),£) + &(£) using Taylor formula, so we have 
6(0) = d^b(0) = and we split the integral (fT5|) with (V6(0),£) and &(£). Since the integrand 
with the (V6(0), £) is smooth and compactly supported in £ (recall that xi = on the support 
of d z x), we can apply stationary phase to get that 



d z x(z)e h 



2^(0 

e h — 

C »• 



£ 



Xi(O<V6(0U> 



7T 



< Ch 2 



uniformly in z. Now set b z (£) = d z x(z)xi(£,)b(0 / ( z ~ which is C 2,a in £ and smooth in z. 
Let 6 £ Cq°([0, 1)) be a cutoff function which is equal to 1 near and set 0h{Q '■= Q{\£\/h), 
then we have by integrating by parts 



(16) 



2«KO~ 
e h b z ^)d^dC 2 



supp(xi) 



supp(xi) 




Using polar coordinates with the fact that b z (0) = 0, it is easy to check that the second 
term in (I16p is bounded uniformy in z by Ch 2 . To deal with the first term, we use 6 2 (0) = 
d^b z (0) = dcb z (0) = and a straightforward computation in polar coordinates shows that the 
first term of (I16p is bounded uniformly in z by Ch 2 \ log(/i)|. We conclude that 

IMIl* <C||7/|| L cc < C^|log h\. 

It is also direct to see that the same estimates holds with a loss of h~ 2 for any derivatives in 
z, z of order less or equal to 2, since they only hit the x( z ) factor, the (z — £) _1 factor or the 
oscillating term e - 2l/ 4>( z )/ h . So we deduce that 

\\ V \\ H2 =0(\\ogh\). 
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and this ends the proof. □ 

We summarize the result of this section with the following 

Lemma 4.3. Let k G N be large and <3? G C k (Mo) be a holomorphic function on Mq which is 
Morse in Mq with a critical point at p G int(Mo). Let a G C k (Mo) be a holomorphic function 
on Mq vanishing to high order at any critical point of other than p. Then there exists 
r% G H 2 (Mq) such that ||ri||£,2 = 0(h) and 

e -*/fc(A + V)e*/ h (a + n) = L2 (fc| log h\). 

4.2. Construction of dp. We have constructed the correction terms ri which solves the 
Schrodinger equation to order h as stated in Lemma l4.3l In this subsection, we will construct 
a holomorphic function ao which annihilates the boundary value of the solution on Tq. In 
particular, we have the following 

Lemma 4.4. There exists a holomorphic function ao G H 2 (Mq) independent of h such that 
e -*/h(A + V)e*' h (a + ri + ha ) = L *{h\ logh\) 

and 

[e* //l (a + n + ha ) + e*/ ft (a + n + /ia )]|r = 0. 

Proof. First, notice that h~ 1 ri\gM = ri2\dM £ H 3 ' 2 (8Mq) is independent of ft. Since $ 
is purely real on Tq and a is purely imaginary on Tq, we see that this Lemma amounts to 
construct a holomorphic function ao G H 2 (Mq) with the boundary condition 

R.e(fi 2 ) + Re(a ) = on Tq. 

To construct ao, it suffices to use (ii) in Corollary 12.21 □ 

4.3. Construction of r<i- The goal of this section is to complete the construction of the 
complex geometric optic solutions by the following proposition: 

Proposition 4.1. There exist solutions to (A + V)u = with boundary condition u\r = 
of the form f 1 1 [> with r\, oq constructed in the previous sections and r 2 satisfying ||t*2 [| x, 2 = 
0(h 3 / 2 |logh|). 

This is a consequence of the following Lemma (which follows from the Carleman estimate 
obtained above): 

Lemma 4.5. IfV& L°°(Mq) and f G L 2 (Mq), then for all h > small enough, there exists 
a solution v G L 2 to the boundary value problem 

e^ h (A g + V)e~^ h v = f, v\ To = 0, 

satisfying the estimate 

\\v\\ L 2<Ch-2\\f\\ L2 . 

Proof. The proof is the same as Proposition 2.2 of [16j . we repeat the argument for the 
convenience of the reader. Define for all h > the real vector space A := {u G Hq(ALq); (A g + 
V)u G L 2 (Mq), d u u |r= 0} equipped with the real scalar product 

(u,w) A := [ e~ 2lp/h (A g u + Vu)(A g w + Vw)dvg. 

J M 
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Observe that since ip is constant along To, d v <p = on Tq. Therefore, we may apply the 
Carleman estimate of Proposition (|3.ip to the weight ip to assert that the space A is a 
Hilbert space equipped with the scalar product above. By using the same estimate, the linear 
functional L : w — > f M e~ ip / h fw dv g on A is continuous and its norm is bounded by h? \ \f\ \ L 2 . 
By Riesz theorem, there is an element u G A such that (., u)a = L and with norm bounded by 
the norm of L. It remains to take v := e- v/h (A g u + Vu) which solves (A g + V)e~ Lp v = e~fl h f 
and which in addition satisfies the desired norm estimate. Furthermore, since 



I e- v/h v(A g + V)wdv g = [ e- v/h fwdv t 

J Mn JMn 



I M J M 

for all w G A, we have by Green's theorem 

^ h vd u wdv g = = [ e- v/h vd u wdv g 
8M Jt 

for all w £ A. This implies v = on Tq. □ 
Proof of Proposition 14.11 We note that 



(A + V)(e* /h (a + n + hao) + e*/ h (a + n + ha ) + e v/h r 2 ) = 

if and only if 

e -v/ft(A + V)e Lp/h r 2 = -e-^ /h (A + V)(e* /h {a + n + ha ) + e^a + n + ha )). 

By Lemma 14.41 the right hand side of the above equation is 0^2 (h\ log Therefore, using 
Lemma 14.51 one can find such r 2 which satisfies 



r 2 



\ L 2 < C^llog h\, r 2 \ ro =0. 



Since the ansatz e^l h {a + n + /iao) + e^l h {a + n + /iao) is arranged to vanish on To, the 
solution 



u = e * /h (a + n + ha ) + e®/ h (a + n + /ia ) + e v/ V 2 
vanishes on Tq as well. □ 



5. Identifying the potential 

We now assume that V\,V 2 G C 1,0! (Mo) are two real valued potentials, with a > 0, such 
that the respective Cauchy data spaces 6^,62 for the operators A g + V\ and A g + V 2 on 
r C 9Mo are equal. Let To = dM$ \ T be the complement of T in <9Mo, and possibly by 
taking V slightly smaller, we may assume that To contains an open set. Let p G M$ be an 
interior point of Mq such that, using Proposition 12, 1\ we can choose a holomorphic Morse 
function $ = + iip on Mo with $ purely real on To, C fc in Mo for some large k G N, with 
a critical point at p. Note that Proposition 12.11 states that we can choose such that none 
of its critical points on the boundary are degenerate and such that critical points do not 
accumulate on the boundary. 

Proposition 5.1. If the Cauchy data spaces agree, i.e. if = CJ> , then V\{p) = V 2 (p). 
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Proof. Let a be a holomorphic function on Mq which is purely imaginary on Tq with a(p) 7^ 
and a(p') = to large order for all other critical point p' of The existence of a is insured 
by Lemma 12.31 Let U\ and u 2 be H 2 solutions on Mq to 

(A g + V 3 ) Uj = 

constructed in Section 2] with $ = (f) + iip for Carleman weight for u\ and — <3? for u%, thus of 
the form 



ux = e® /h (a + ha + r 1 ) + e*/ h (a + ha + n) + e (p/h r 2 



u 2 = e~^ h (a + hbo + si) + e-*A(a + /i6 + si) + e~^ h s 2 

and with boundary value Uj\gM = fj, where fj vanishes on To- We can write by Green 
formula 



/ Ul(Vl - V 2 )U2dVg = - I (AgUl.U 2 - Ul.AgU 2 )dVg 

J M J M 

= - (d u ui.f 2 - fi.d v u 2 )dv g . 

JdMn 



Id M 

Since the Cauchy data for A g + V\ agrees on V with that of A g + V 2 , there exists a solution 
v of the boundary value problem 

(A g + V 2 )v = 0, v\ 9Mo = fx, 

satisfying d u v = d v u\ on T. Since fj = on Fq, this implies that 

(17) 

/ ui(Vx - V 2 )u 2 dv g = - / (A g ui.u 2 - ui.A g u 2 )dvg = - / (d v ui.f 2 - f\.d v u 2 )dv g 

J M J Mo J 8 Mo 

= - / (d p v.f 2 - v.d v u 2 )d\ g = — (A g v.u 2 - v.A g u 2 )dv g = 
JdM Jm 

since A g + V 2 annihilates both v and u 2 . We substitute in the full expansion for u\ and u 2 
and, setting V := V\ — V 2 , and using the estimates in Lemmas 14.21 14. l l and 14.41 and, we obtain 

(18) = I x + h + o(h), 
where 

(19) h= [ V(a 2 + a 2 )dv g + 2Re [ e 2i ^ /h V\a\ 2 dv g , 

J M Q J Mo 

(20) 1 2 = 2hRe [ aV(e 2i ^ h (^ + b ) + e- 2i ^ h (a + + &o + «o + ^^)dv 9 . 

Jm \ n n n / 

Remark 5.1. W^e observe from the last identity in Lemma \4-S\ that r\/h in the expression I 2 
can be replaced by the term r\ 2 satisfying 2iri 2 dtp = b up to an error which can go in the o(h) 
in ()18p . and similarly for the term S\/h which can be replaced by a term S12 independent of 
h. 

We will apply the stationary phase to these two terms in the following two Lemmas. 
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Lemma 5.2. The following estimates holds true 

I 2 = 2hRe( / Va(b + a + r 12 + S12 + S12 + rn)dvg) + o(h). 
V J M ' 

where ryi, syi, f\i ands\2 are independent of h. 

Proof. We start by the following 

Lemma 5.3. Let f G L 1 (Mo), then as h — > 

e ^/h fdVg = o(1) 

Proof. Since C k (Mo) is dense in L 1 (Mo) for all k G N, it suffices to prove the Lemma for 
/ G C k (Mo). Let e > be small, and choose cut off function \ which is identically equal to 
1 on the boundary such that 

/ Xl/|dv 9 <e. 
J Mo 

Then, splitting the integral and using stationary phase for the 1 — % term, we obtain 



Mo 



< 



(1 - x)e 2 ^ /h fdv 



Mo 



+ 



Mo 



<e + O e (h) 



which concludes the proof by taking h small enough depending on e. 

The proof of Lemma 15.21 is a direct consequence of Lemma 15.31 and Remark 15.1 



□ 
□ 



The second Lemma will be proved in the end of this section. 
Lemma 5.4. The following estimate holds true 

I 1= [ V{a 2 + a 2 )dv g + hC p V(p)\a(p)\ 2 Re{e 2i ^ p V h ) + o(h) 
JMo 

with C p 7^ and independent of h. 

With these two Lemmas, we can write f)18|) as 

0= f V(a 2 + a 2 )dv g + 0(h) 
JMo 

and thus we can conclude that 

0= f V(a 2 + a 2 )dv g . 
JMo 

Therefore, (fl8|) becomes 

= C p V(p)\a(p)\ 2 Re(e 2l ^/ h ) + 2Re( / Va(b + Sl + a + n)dv g ) + o(l). 

v J Mo ' 

Since tp(p) 7^ we may choose a sequence of hj -> such that Re(e 2i ^PV h i) = 1 and another 

sequence hj — > such that Re(e 2 *^^ ///i J ') = — 1 for all j. Adding the expansion with h = hj 

and h = hj, we deduce that 

= 2C p V(p)\a(p)\ 2 + o(l) 
as j — ► 00, and since C p 7^ 0, a(p) 7^ 0, we conclude that V(p) = 0. The set of p G Mq for 
which we can conclude this is dense in Mq by Proposition 12.11 Therefore we can conclude 
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that V{p) = for all p E Mq. □ 



We now prove Lemma 15.41 
Proof of Lemma 15.41 Let x De a smooth cutoff function on Mq which is identically 1 
everywhere except outside a small ball containing p and no other critical point of and 
X = near p. We split the oscillatory integral in two parts: 

( e 2#M + e - 2 ^ /h )V\a\ 2 dv g = f X {e 2ilp/h + e- 2i ^ h )V\a\ 2 dv g 
Mo J M 

+ [ (1 - x)(e 2#/ft + e- 2i ^ h )V\a\ 2 dv g 

J Mo 

The phase ip has nondegenerate critical points, therefore, a standard application of the sta- 
tionary phase at p gives 

/ (1 - x)(e 2iij/h + e- 2i ^ /h )V{p)\a\ 2 dv g = hC p \a(p)\ 2 V(p)Re{e 2i ^ p)/h ) + o(h) 

Jm 

where C p is a non-zero number which depends on the Hessian of ip at the point p. Define the 
potential V{-) := V(-) - V(p) G C 1,Q (M ), then we show that 

(21) / (1 - X )(e 2i ' p/h +e- 2ii,/h )V\a\ 2 dv g = o(h). 

J Mo 

Indeed, first by integration by parts and using A g tjj = one has 

(l-x)(e 2 ^ /h + e- 2l ^ h )V\a\ 2 dv g =^ [ (d(e 2l ^ h - e~ 2 ^ h ), d^V^j^dv, 



h f (e 2#A _ e -**/h ){d ( 0- -x]\f V ), #)dv fl 

M 



2i 



but we can see that (d((l - x)\a\ 2 V/\dip\ 2 ),dip) G L 1 (M ): this follows directly from the 
fact that V is in the Holder space C 1,cx {Mq) and V(p) = 0, and from the non degeneracy 
of Hess(V0- ^ then suffice to use Lemma 15.31 to conclude that (|21h holds. Using similar 
argument, we now show that 

f x{e 2ilp/h + e- 2i ^ /h )V\a\ 2 dv g = o(h). 
J Mq 

Indeed, since a vanishes to large order at all boundary critical points of ip, we may write 



x(e W + e -2#M )F | a |2 dv " / {d{e 2i^/h_ e -2i^ d ^ v Xp_ dVg 

M 2^ J Mq 



2% 

For the interior integral we use Lemma 15.31 to conclude that 

^ f rSiiblh „-2iiblh\ firXW 2 
'Mo 



h f _ e -**/A )di Yg ( V ^% W 

2« Jm v l#| 2 

/ {e ^/h_ e -^/h )v ^_ dulpdv ^ 

dM 



2i 



I 1 2 
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and for the boundary integral, we write 8Mq = Tq U T and observe that on To, tp = so 
^ e 2iip/h _ e -2itj}/h^ _ whiig on y we have V = from the boundary determinacy proved in 
Proposition 17.11 of the Appendix. Therefore 

f X{e m,h + e- 2i ^ h )V\a\ 2 dv g = o(h) 

and the proof is complete. □ 



6. Inverse scattering 

We first obtain, as a trivial consequence of Theorem 11.11 a result about inverse scattering 
for asymptotically hyperbolic surface (AH in short). Recall that an AH surface is an open 
complete Riemannian surface (X, g) such that X is the interior of a smooth compact surface 
with boundary X, and for any smooth boundary defining function x of dX, g := x 2 g extends 
as a smooth metric to X, with curvature tending to —1 at dX. If V G C°°(X) and V = 
0(x 2 ), then we can define a scattering map as follows (see for instance |17l I12j or [13j): first 
the L 2 kernel kev L 2(A g + V) is a finite dimensional subspace of xC°°(X) and in one-to-one 
correspondence with E := {(d x ^)\gx'i V> £ ker^A^ + V)} where d x := V% is the normal 
vector field to dX for g, then for / G C°°(dX), there exists a function u £ C°°(X), unique 
modulo ker^2(A s + V), such that (A s + V)u = and u\qx = f ■ Then one can see that the 
scattering map S : C°°(dX) — > C°°(dX)/E is defined by 8/ := d x u\ dx . We thus obtain 

Corollary 6.1. Let (X,g) be an asymptotically hyperbolic manifold and let V\, V2 £ x 2 C°°(X) 
be two potentials and V C dX an open subset of the conformal boundary. Assume that 

{9xu\ d x;u G ker i2 (A 9 + Vi)} = {d x u\ 9x ;u G ker L2 (A 9 + V 2 )} 

and let §>j be the scattering map for the operator A g + Vj for j = 1, 2. 7/§i/ = §2/ on T for 
all f G C£°(r) then Vi = V 2 . 

Proof. Let x be a smooth boundary defining function of dX, and let g = x 2 g be the 
compactified metric and define Vj := Vj/x 2 G C°°(X). By conformal invariance of the 
Laplacian in dimension 2, one has 

A g + V j = x 2 {A- 9 + V j ) 

and so if ker L 2(A 5 + V\) = ker L 2(A 9 + V2) and §i = §2 on T, then the Cauchy data spaces 
Cf for the operator Ag+Vj are the same. Then it suffices to apply the result in Theorem ll.il □ 

Next we consider the asymptotically Euclidean scattering at frequency. An asymptotically 
Euclidean surface is a non-compact Riemann surface (X,g), which compactifies into X and 
such that the metric in a collar (0, e) x x dX near the boundary is of the form 

dx 2 h{x) 



where h{x) is a smooth one-parameter family of metrics on dX with /i(0) = d0 2 sl is the metric 
with length 2ir on each copy of S 1 that forms the connected components of dX. Notice that 
using the coordinates r := 1/x, g is asymptotic to dr 2 + r 2 d6 2 sl near r — ► 00. A particular 
case is given by the surfaces with Euclidean ends, ie. ends isometric to M 2 \ B(0,R) where 
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-6(0, R) = {z E R 2 ; \ z\ > R}. Note that g is conformal to an asymptotically cylindrical metric, 
or 'b-metric' in the sense of Melrose |24j . 

g b :=x g = —j- + /i(x) 

and the Laplacian satisfies A g = x 2 A gb . Each end of X is of the form (0, e) x x and the 
operator A gb has the expression in the ends 

A 9t = -(xdx) 2 + A d x + xP(x,8;xd x ,d e ) 

for some smooth differential operator P(x,9;xd x ,dg) in the vector fields xd x ,dg down to 
x = 0. Let us define V b := x~ 2 V, which is compactly supported and 

H 2m :={uG L 2 (X, dvol Si) ); A™«£ L 2 (X, dvolj}, m E N . 

We also define the following spaces for a E R 

F a := kei xaH 2(A gb + V b ). 

Since the eigenvalues of A s i are {j ,j E No}, the relative Index theorem of Melrose \24\ 
Section 6.2] shows that A gb + 14 is Fredholm from x a H 2 to x a H® if a Z. Moreover, 
from subsection 2.2.4 of [23], we have that any solution of (A gb + V b )u = in x a H 2 has an 
asymptotic expansion of the form 

u~ '^2xi (log x) £ Uj ^(6), as x — > 

for some sequence of non negative integers and some smooth function Uj ; £ on S . In 
particular, it is easy to check that kei L 2^ X) dvoig)(^g + V) = Fi+t for e E (0, 1). 

Theorem 6.2. Let (X,g) be an asymptotically Euclidean surface and Vi, V% be two compactly 
supported smooth potentials and x be a boundary defining function. Let e E (0, 1) and assume 
that for any j E Z and any function ip E ker x j- eH 2 (A g + V\) there is a tp E ker a ,j- e ^2 (A g + V2) 
suc/i that ip — ip = 0(x°°), and conversely. Then V\ = V%. 

Proof. The idea is to reduce the problem to the compact case. First we notice that by 
unique continuation, ip = ip where V\ = V% = 0. Now it remains to prove that, if R v denote 
the restriction of smooth functions on X to {x > rj} and V is a smooth compactly supported 
potential in {x > n}, then the set U°? =0 i?, ? (F__j„ e ) is dense in the set Ny of H 2 ({x > n}) 
solutions of (A g + V)u = 0. The proof is well known for positive frequency scattering (see 
for instance Lemma 3.2 in [25]), here it is very similar so we do not give much details. The 
main argument is to show that it converges in L 2 sense and then use elliptic regularity; the 
L 2 convergence can be shows as follows: let / E Ny such that 

f fipdvolg = 0, E Uj°- F-j- e , 

then we want to show that / = 0. By Proposition 5.64 in [24] , there exists k E N and a 
generalized right inverse G b for P b = A gb + V& (here, as before, x 2 V b = V) in x~ k ~ e H 2 , such 
that PbGb = Id. This holds in x~ k ~ e H 2 for k large enough since the cokernel of P& on this 
space becomes for k large. Let u> = G b f so that (A gb + Vh)uj = f, and in particular this 
function is in {x < rf\. The asymptotic behaviour of the integral kernel Gb(z,z') of G b as 
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z — ► oo is given in Proposition 5.64 of [21] uniformly in z' £ {x > rj}, we have for all J £ N 
and using the radial coordinates (x, 6) £ (0, e) x S 1 for z in the ends 

J ij 

G b (z,z') = ]T ^2^(logx)%(e,z')+o(x J ) 

j=-k i=0 

for some functions ipj£ £ x k ~i~ e H 2 and some sequence (£j)j of non-negative integers. But 
the fact that (A 9b + Vb)Gb(z, z') = 5(z — z') as distributions implies directly that (A gb + 
Vb)i/jj(9,.) = 0. Using our assumption on /, we deduce that J x ipj(0, z')f(z')dvol gb = for 
all j £ No and so the function u vanish faster than all power of x at infinity. Then by unique 
continuation, we deduce that uj = in {x < e}. Since now u £ H 2 , its Cauchy data at x = rj 
are and A gb + Vj, is self adjoint for the measure dvol gb , we can use the Green formula to 
obtain 

f \f\ 2 dvol gb = f u(A gb + V b )fdvol gb = 0. 

J X>7] J X>T] 

The H 2 density is easy using elliptic regularity. □ 



7. Appendix : boundary determination 

In this appendix, we give a short proof of the fact that the partial Cauchy data space on 
r C dM determines the potential on V when the potential is in C°' a (M) for some a £ (0, 1). 
This result is shown for the conductivity problem on a domain of R n in [19] under the 
assumption that the conductivity has roughly n/2-derivatives, it is also shown in [29J for 
continuous potentials on a smooth domain of W 1 by using a limiting argument from the smooth 
case (which they analyze using micolocal analysis near the boundary). Alessandrini [lj also 
proved such a result under the assumption that the domain is Lipschitz and the coefficients 
of the operator are in W 1,p for p > n, while Brown [H] studied the case of Lipschitz domains 
with a continuous conductivity. Since the result in our setting is not explicitly written down, 
but certainly known from specialists, we provide a short proof without too many details, and 
using the approach of [5] . We shall prove 

Proposition 7.1. Let T C 8Mq be an non-empty open subset of the boundary. If V\,Vi £ 
C°' a (M) for some a > and their associated Cauchy data spaces C5\ C?> defined in (fT]) are 
equal, then V\\r = V-^r- 

The key to proving this proposition is the existence of solutions to (A g + Vi)u = which 
concentrate near a point p £ T. First we need a solvability result for the equation (A g + Vi)u = 
f, which is an easy consequence of the Carleman estimate of Proposition 13.11 and follows the 
method of Salo-Tzou [27} Section 6]. If we fix h > small and take ip = 1 in the Carleman 
estimate of Proposition 13.11 we obtain easily that there is a constant C such that for all 
functions in H 2 (Mq) satisfying u\qm = 

(22) Hik + R«||f 2(ro) < a(||(A + ^H|| 2 + ||M| 2(r) ) 

As a consequence, we deduce the following solvability result: let 

25 := {w £ H 2 (M ) n f#(M ) | d„w\ T = 0} 
be the closed subspace of H 2 (Mq) under the H 2 norm and let 23* be its dual space then 
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Corollary 7.1. Let i = 1,2, i/ien /or a// / G L 2 (Mq) there exists u G H 2 (Mq) solving the 
equation 

(Ag + ^)« = Z 

witft boundary condition u\r = 0, and \\u\\jj2 < C||/||b*- 

Proof. Set yi := {w G i?o(Mo) | (A g + V^)w G L 2 ,e?,/u;|r = 0} equipped with the inner 
product 



:= / (A g + Vi)v(A g + Vi)wdv g . 

J Mo 



(V,W)A 

'M 

Thanks to ([22]) . A is a Hilbert space and A = S. For each / G L 2 (Mo), let us define the 
linear functional on 13 



Lf : w i— > / io/ dv g . 

./Mo 



By fl22J) , we have that for all w G S 

\L f (w)\ < H/b'Ikb < 11/b'lklU- 

Therefore, by Riesz Theorem, there exists Vf € A. such that 



/ (A s + y)u7(Ag + Vi)u;dv 9 = / ^/d\ 

J M J M 



for all it; G 25. Furthermore, ||(A ff + T^u/H^a < ||/||b*- Setting u := (A g + Vi)vJ we have that 
(A g + Vi)u = / and ||w||i2 < To obtain the boundary condition for u, observe that 

since 



/ u(A g + Vi)wdvg = / fwdvg 

JMn J Mo 



I Mo J Mo 

for all w G 23, by Green's theorem 



ud v w dv 9 = = / ud u w dv 5 



'Mo J To 

for all to G 23. This implies it = on Tq. □ 



Clearly, it suffices to assume that T is a small piece of the boundary which is contained 
in a single coordinate chart with complex coordinates z = x + iy where \z\ < 1, Im(^) > 
and the boundary is given by {y = 0}. Moreover the metric is of the form e 2p \dz\ 2 for some 
smooth function p. Let p G T and possibly by translating the coordinates, we can assume 
that p = {z = 0}. Let r\ G C co (Mq) be a cutoff function supported in a small neighbourhood 
of p. For h > small, we define the the smooth function Vh G C°°(Mq) supported near p via 
the coordinate chart Z = (x, y) G M 2 by 

(23) v h (Z) := n(Z/Vh)e$°" z 

where a := (i, —1) G C 2 is chosen such that a ■ a = 0. We thus get (d 2 + d 2 )e a ' z = and thus 
A g e a ' z = by conformal covariance of the Laplacian. Therefore, we have in local coordinates 

(24) A g v h (Z) = l e^ z (A a r,)(^=) + J- e ^ z (dr,A),a.dZ) g . 
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Lemma 7.2. IfV £ C 0,a (M ) for q > 2, then there exists a solution Uh G H 2 to (A g + V)u = 
of the form 

Uh = v h + R h , 

with Vh defined in ([23]) and WR^W^ < Ch^^ , satisfying swpp(Rh\dM ) C T. 

Proof of Lemma 17.21 We need to find Rh satisfying HU^Hi^ < C/i 5//4 and solving 

(A g + V)R h = -(A g + V)v h =: M h . 

Thanks to Corollary 17,11 it suffices to show that H-M/JIb* < C/i 5//4 . Thus, let w & 23, then we 
have by flSD 

/ wM h dv g = I x + I 2 + h 

J Ma 



where 



L 



wViV h e 2 f>dZ, h:=~ [ we^ z X i{^f)e 2p dZ, 

; n J\z\<vK vh 

[ WX2 A)e*'* z dZ 
J\z\<Vh Vh 



h 3 / 2 J\z\<Vh 

and xi = A 9 r/, Xi = id x n — d y r\. In the above equation the third term ^3 has the worst growth 
when h — ► 0. We will analyze its behavior and the preceding terms can be treated in similar 
fashion. One has 

h = ~ h 1 ' 2 f WX2 (^)(d x - id y ) 2 e^ z dZ 
J\z\<Vh Vh 

= ~ h 1 ' 2 f (d x - id y ) 2 {wx2(^=))e^ x dZ. 
J\z\<VJi Vh 

Notice that the boundary term in the integration by parts vanishes because w £ Hq and 
9 v w\dM vanishes on the support of rj. The term (d x — id y ) 2 ( W X2(-^)) has derivatives hitting 
both x2(-^ h =) and w. The worst growth in h would occur when both derivatives hit X2(^0 i n 

which case a h^ 1 factor would come out. Combined with the h 1 / 2 term in front of the integral 
this gives a total of a h~ l l 2 in front. By this observation we have improved the growth from 
/i~ 3 / 2 to ft, -1 / 2 . Repeating this line of argument and using Cauchy-Schwarz inequality, we can 
sec that \h\ < Ch 5 / 4 \\w\\ H 2 (an elementary computation shows that functions of the form 
X {Z/Vh)e^ z have L 2 norm bounded by Ch 3 ' 4 ). Therefore, H^^^e^b' < C/1 5 / 4 
and we are done. □ 



Proof of Proposition [7TT1 It suffices to plug the solutions ul,ul from Proposition 17.21 into 
the boundary integral identity (|17fl . A simple calculation using the fact that V\ — V 2 is in 
C7°' 7 (M ) yields that 

= / utiVi - V 2 )u 2 h dv 9 = Ch^iVxip) - V 2 {p)) + o(h 3 / 2 ) 
J M 

and we are done. □ 
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